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Major  Research  Findings 

The  original  proposal  contained  four  parts,  and  we  have  divided  the  final  report  in  the  same  way.  In 
each  part,  we  briefly  describe  our  major  findings. 

1.  Flow  directed  iterative  methods 

It  was  proposed  to  develop  iterative  methods  for  non-symmetric  matrices  of  the  type  arising  from 
convection-diffusion  problems. 

Our  results  in  this  area  are  contained  in  [1,4,5].  [1]  deals  with  an  iterative  method  for  a  particular  prob¬ 
lem,  and  illustrates  the  usefulness  of  shaping  the  iterative  method  to  the  problem  at  hand.  Our  major  result 
in  this  area  is  contained  in  [4],  in  which  it  is  shown  that  for  convection  diffusion  problems,  an  appropriate 
choice  of  iteration  parameters  gives  work  estimates  for  the  solution  of  the  problem  that  are  comparable  to 
work  estimates  for  the  model  elliptic  self  adjoint  problem;  that  is,  optimally  over-relaxed  SOR  applied  to 
the  model  5  point  formula  for  the  Laplace  operator.  In  [5]  we  have  given  a  survey  of  some  of  our  results 
and  conjectures  in  this  area.  We  are  planning  some  numerical  experiments  to  test  these  results  in  the  near 
future. 

2.  Numerical  solution  of  parabolic  problems  with  interior  layers 

The  proposed  work  consisted  of  studying  the  asymptotic  structure  of  internal  layers  and  moving  fronts. 
Our  results  in  this  area  have  not  achieved  final  form,  so  we  have  not  published  anything  on  this. 
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3.  Viscous  compressible  flow 

The  proposed  work  in  this  area  is  to  establish  the  existence  of  discontinuous  viscous  compressible  flows, 
and  to  determine  important  properties  of  these  flows. 

Our  work  in  this  area  is  contained  in  [3,6].  Each  of  these  papers  contains  major  findings.  In  [3]  we 
have  established  the  existence  of  discontinuous  solutions  of  a  nonlinear  system  of  partial  differential  equations 
with  an  elliptic-hyperbolic  character.  The  system  is  related  to,  but  simpler  than  the  Navier-Stokes  equations. 
The  analysis  depends  on  the  Ccireful  choice  of  Banach  space;  the  norm  in  this  space  must  be  strong  enough 
to  fulfill  demands  set  by  the  hyperbolic  part  of  tbe  system,  and  weak  enough  to  fulfill  demands  set  by  the 
elliptic  part  of  the  system.  This  work  constitutes  the  Ph.D  thesis  of  Senhuei  Chen,  who  was  partly  supported 
by  the  grant. 

In  [6]  we  give  a  closer  study  of  the  linearized  viscous  compressible  flow  equations,  linearized  at  a  uniform 
ambient  flow.  Using  the  Fourier  transform  and  Mathematical  the  most  important  properties  of  a  discontin¬ 
uous  solution  are  found.  These  properties  include  an  rln  r  type  of  singularity  in  the  velocity  at  the  entrance 
and  exit  points  of  the  curve  of  discontinuity.  These  results  are  being  rewritten  for  publication  in  a  refer¬ 
eed  journal.  We  have  continued  this  line  of  investigation,  cuid  have  formulated  a  linearization  of  the  flow 
equations  around  an  ambient  discontinuous  solution  of  tbe  nonlinear  system.  Using  this  we  have  established 
the  linecirized  stability  of  putative  discontinuous  viscous  compressible  flows  that  are  close  to  uniform.  These 
results  are  also  being  written  up  for  publication. 

4-  Corners  in  singular  perturbation  problems 

The  proposed  work  here  is  to  study  the  interaction  of  boundary  layers  and  corner  singularities  in 
singularly  perturbed  elliptic  boundary  value  problems  in  polygonal  domains. 

[2]  contains  our  published  work  in  this  area.  We  have  actively  continued  this  work,  and  obtained 
major  new  results  that  will  shortly  be  available.  (A  copy  will  be  sent  to  the  ARO  when  it  is  submitted  for 
publication.)  The  new  work  considers  the  singularly  perturbed  elliptic  boundary  value  problem 

Lu  :=  —e^Au  +  qu  =  f  in  Cl, 
u  =  g  onT  =  dCi, 

where  fi  C  is  the  interior  of  a  polygon  F.  ([2]  considers  the  case  of  a  rectangle.)  The  study  involves  the 
interaction  of  two  distinct  phenomena,  corner  singularities  and  boundary  layers.  We  obtain  an  asymptotic 
expansion  of  the  solution  that  displays  both  the  corner  singularities  occuring  near  the  vertices  of  F  and  the 
boundary  layer  behavior  that  is  present  when  e  is  small.  The  asymptotic  expansion  can  be  differentiated  to 
provide  an  asymptotic  expansion  for  the  derivatives  of  u.  As  a  consequence  of  the  asymptotic  expansion, 
bounds  for  tbe  derivatives  of  the  solution  are  derived  that  display  both  the  dependence  on  e  and  the  depen¬ 
dence  on  the  distance  to  a  corner  or  side  of  F.  These  bounds  show  that  while  a  small  value  of  s  produces 
boundary  layers  along  the  sides  of  F  that  decay  exponentially  as  one  moves  into  Q,  the  small  value  of  e  also 
mitigates  the  effect  of  the  corner  singularities.  That  is,  the  corner  singularity  is  multiplied  by  an  attenuation 
factor  of  the  form  where  d„  is  the  distance  to  the  corner.  This  means,  for  example,  that  the  pollution 

effect  of  the  corner  singularity  on  the  error  in  a  finite  element  computation  will  be  mitigated  by  a  small  value 
of  s. 
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